Time evolution of correlation functions in Non-equilibrium Field Theories 



o 

Q 

oo 



> 

^ . 

Oh 
I 

Oh 



X 



Luis M. A. Bettencourt^'^ and Christof Wetterich^ 
^ Institut fur Theoretische Physik, University of Heidelberg, Philosophenweg 16, 69120 Heidelberg, Germany 
^T6/T11, Theoretical Division, MS B288, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

(February 1, 2008) 

We investigate the non-equilibrium properties of an Ai'-component scalar field theory. The time 
evolution of the correlation functions for an arbitrary ensemble of initial conditions is described 
by an exact functional differential equation. In leading order in the 1/N expansion the system 
can be understood in terms of infinitely many conserved quantities. They forbid the approach 
to the canonical thermal distribution. Beyond leading order only energy conservation is apparent 
generically. Nevertheless, we find a large manifold of stationary distributions both for classical 
and quantum fields. They are the fixed points of the evolution equation. For small deviations of 
the correlation functions from a large range of fixed points we observe stable oscillations. These 
results raise the question of if and in what sense the particular fixed point corresponding to thermal 
equilibrium dominates the large time behavior of the system. 
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The understanding of the time evolution for a quantum 
field theory with initial conditions away from thermal 
equilibrium is needed for a large variety of phenomena, 
ranging from inflationary cosmology in the early universe 
to laboratory experiments of defect formation. First in- 
sights have been gained by perturbative methods and, in 
particular, by an investigation of the leading order in the 
1 /N expansion, with N the number of components of the 
field ||l| ^ . This approximation has been applied to a va- 
riety of problems like decoherence of quantum states , 
heating after cosmological inflation disoriented chi- 
ral condensates in heavy ion collisions [^, decoherence 
and photon pair creation from a homogeneous electric 
field H]. In most of this work a particle number basis 
was used in order to describe the evolution of each field 
mode. From the mode evolution the two point function 
can in turn be constructed. 

It is a severe limitation, however, that the leading 1/iV 
approximation neglects scattering and cannot describe 
the approach of a system to a thermal equilibrium state. 
A systematic extension of the particle number formal- 
ism to finite TV seems a difficult task. We pursue here 
an alternative strategy and investigate directly the time 
evolution of the correlation functions of the system. One 
starts, at t = 0, with an ensemble of initial values for the 
field variables and their time derivatives. This specifies 
the initial correlation functions. At any later time t the 
correlation functions are, in principle, computable from 
the microscopic field equations. In practice, we follow 
this time evolution by means of an exact functional dif- 
ferential evolution equation for the generating functional 
of the correlation functions Approximate solutions 
can be obtained by truncation. 

Our main tool is the time dependent effective action 
r[(/), 7r;t] 1^. For field theories F is a functional of the 
fields (l){x) and ■k{x) — dt(j){x) and generates the one par- 
ticle irreducible equal time correlation functions. Its time 



evolution obeys 



(1) 



where the first operator, Cc\, generates the classical evo- 
lution, whereas £q determines the dynamics of the quan- 
tum effects. We will deal with the specific example of 
a 0{N) symmetric field theory for real scalars Xa{x), 
a — l,...,iV, obeying the microscopic equations of mo- 
tion of a 0'*-model 



dtXa = did.Xa - m'^Xa - 7:XbXbXa 



(2) 



The arguments of F are the mean values of the field (jja — 
(Xa) and its time derivative tTq = (xa) (in the presence 
of sources). From iGlfTll, one infers 
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with '0 G {4>^ Tr} and repeated index summations implied. 
The exact field dependent propagator G, interpreted as 
a matrix with multi- indices (V', i,x), is the inverse of the 
second functional derivative of F, 
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We emphasize that the minimum of T[(j), tt; t] describes 
the time evolution of a cosmological scalar field like 
the inflaton. The true dynamical equation (|l|) for non- 
equilibrium situations differs, in general, from the varia- 
tion of some equilibrium effective action or the classical 
action. The latter are at best an approximation to Eq. (|l|) 
in certain limiting situations not too far from equilibrium. 
An understanding of the structure of the dynamical equa- 
tions following from (|^) is of crucial importance for an 
understanding of cosmological models like inflation. 

We assume an 0(A^)-symmetric initial distribution for 
which all terms with odd powers of the fields in T van- 
ish. This property remains conserved by the evolution. 
Similarly, we consider a translational and rotational in- 
variant situation. The most general form of F consistent 
with these symmetries can be expanded in terms of the 
Fourier components of the fields 

+2C{q).M)U<l)\ + 8 J (2^)1. (2^)1. (2.)^ 
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Here A, B and C are functions of (7^ and the quartic cou- 
plings M, u, etc., depend on the invariants that can be con- 
structed from the three independent momenta gi, 52, 93. 

The evolution equation can also be expanded in pow- 
ers of the fields and the exact equation for the inverse 
propagator reads 



dtA{q) = 2ojlC{q), 

dtB{q) = -2C{q) - 2^gB{q), 

dtC{q) = LulB{q) - A{q) - j,C{q). 



(7) 



Here, the time dependent frequency, ujq, and the coeffi- 
cient 7g are: 
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and G{q) the exact 



propagator for vanishing fields 
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We suppose that all momentum integrals are ultraviolet 
regularized by a cutoff A, {q^ < A^). For the simple 
system of N coupled anharmonic oscillators {D — 0) the 
momentum labels and integrals should be omitted. From 
(0) one obtains the fiow equations for G{q), c{q) and 
"(9), 



dta^iq) = -2j,a^{q), 
dtG{q) = -2c{q)G{q), 

dtc{q) =ujq - 
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For large N, with A, u, v, w, s, y, z scaling like l/A'^, one 
may expand in XjN as a small parameter. We observe 
that 7g ~ 1/A^ and iOq becomes independent of m, w, . . . , 
in leading order. In the limit N ^ 00 the system of 
evolution equations for the inverse two-point functions 
A, B and C is closed. 

For a solution of Eq. (|^) at finite iV, however, one needs 
the time dependent four point functions. Their evolution 
equations involve, in turn the six point functions. The ex- 
act system is not closed and for any practical purpose we 
have to proceed to some approximation. In the present 
note we simply omit all contributions from 1 PI six point 
vertices. Furthermore for D > 1, we neglect the mo- 
mentum dependence of the quartic couplings. With this 
truncation the evolution equations for the quartic cou- 
plings are given by 
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dtw = 3oj'^y -3v~ 2^w - SXh^B^c 
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For D > 1 there remains at this stage an ambiguity at 
which momentum the couphngs u, v, etc. and therefore 
the barred quantities w^, 7, c, B, should be evaluated. (For 
D — the bars should be omitted.) One may choose 
some appropriate momentum average to be discussed be- 
low. From the truncated evolution equation of the quar- 
tic couplings we have omitted some terms ~ 1 /N'^ which 
involve two-loop integrals. The neglected six point func- 
tions can also be treated consistently as being ^ 1/N^. 
On the other hand, the difference u{q,p,p') — m(0,0,0) 
appears in the contributions ~ 1/7V. Nevertheless, the 
effects of this momentum dependence are partly aver- 
aged out by the momentum integrations. Thus our set 
of evolution equations (|l^) contains for D = all contri- 
butions in the next to leading order This extends 
to arbitrary D only if the momentum dependence of the 
quartic couplings can be neglected. 

The microscopic field equations (||) conserve the total 
energy E for every choice of initial state 



E= I d^'n 



^dtXadtXa + ^diXadiXa + \^rr?XaXa 
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must be conserved by the exact evolution equations. We 
have verified that the truncated evolution equations (|^), 
( p^ ) conserve e for D = Q and arbitrary N . For D > \ 
this holds automatically in leading order in the \/N ex- 
pansion. For finite N our truncation conserves energy 
only provided the averages uP', B, c, 7, are chosen appro- 
priately such that they obey the relation 



dte^ 



N{N + 2] 



X 



d^q d^p d^p' 



■G{q)G{p) 



(27r)^ (27r)^ (27r)^ 

xG{p')G{p + p' + 9)| [(c^' - + (7 - 7,) c(9)] 

[v — 2wc{p) + Zyc{p)c{p') ~ 4:zc{p)c{p')c{p + p' + 9)] 



+4A^ (c-c(9)) l-uSo + -{Si+cip)S2) 



{N + 2)w ~ {N -l)s 
2 {N + 8) 

6w+{N -l)s 



[S2 + c{p)c{p')So] 



(iV + 8) ""^P^^' 2 ^""^P^^^ + c(p)c(p')5i] 



^,S2--h^B^ (c--c(p))(c-c(p')) 



(17) 



For practical purposes it may be sufhcient that this re- 
lation holds in a time averaged sense. Whereas ei is 
separately conserved in the leading 1/iV approximation 
the flow of energy between ei and €2 can be attributed 
to scattering. 

The fixed points of the system Eqs. (|^), ( p^ corre- 
spond to stationary probability distributions. They obey. 



for > 0, 



A4q)=uj^gB4q), 
=0, = 0, 

_ AS* -I- (wV2) 



a(9) =0, 



1 -l- XBsfSo 
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2[jj^z* + j^^g XB^Sqs^ — B^ 
^* ^ 1 I jV+2 ■> • 

We observe a large manifold of fixed points since Eqs. (|l^ ) 
have solutions for arbitrary B^{q), and ! This prop- 
erty seems not to be an artifact of the truncation. It 
persists for momentum dependent four point functions. 
The fixed point manifold becomes even larger if one in- 
cludes the six point functions and seems to characterize 
any higher (finite) polynomial truncation. Classical ther- 
modynamic equilibrium, at temperature T (for h = 0) 
corresponds to the particular point in this manifold: 
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Beq(g) =P= Ceq(g) = 0, 



(19) 



<Jcq 



>Jcq 



Ucq = A/3 (1 ~ SoUcq) . 

Here the fixed point conditions for ^eq and Meq con- 
stitute classical Schwinger- Dyson equations (correspond- 
ing to a Euclidean D-dimensional quantum field the- 
ory at zero temperature) for the effective propagator 
(mass and kinetic terms) and the effective quartic cou- 
pling, respectively. In quantum field theory (h — 1) the 
thermal equilibrium fixed point does not correspond to 
Wcq = Scq — Zcq = 0, auymorc. The Schwinger-Dyson 
equations have now to be evaluated in a D+1 dimensional 
quantum field theory with periodic boundary conditions 
in 'Euclidean time'. The difference between the quan- 
tum and classical fixed points is due to the nonzero Mat- 
subara frequencies. We observe, nevertheless, a rather 
similar structure of the classical and quantum dynami- 
cal systems. In particular, we find in the quantum sys- 
tem a fixed point, which differs from the classical equi- 
librium fixed point only by = Xfi^ / (AH)^) , instead of 
^^oq = 0. For this fixed point, all correlation functions in- 
volving only (f> are the same as for classical equilibrium. 
Inversely, for every fixed point of the quantum system 
there is a corresponding one for the classical system, dif- 
fering only by a shift in z^. Up to this shift, the quantum 
equilibrium could be reached by the classical system ! 

The central issue of this work is the question of whether 
and how the fixed points are approached for t —^ oo. We 
start with the leading behavior for iV — > oo, where 



with 



and 



79 = 0, 



XN 



i2n) 



d'^V C{p) _ 



(20) 



(21) 



(22) 



In this limit the Fourier components of the determinant 
{q) become an infinite set of conserved quantities ! (cf. 
Eqs. ([l^), (p^).) Obviously, the system can reach ther- 
mal equilibrium only for a very particular set of initial 
conditions, namely if a^(g) = G^^{q), for every q. 

Nevertheless, for large times and D > 1, Am^ (or w^) 
tends to a constant value as can be seen in Fig. |l|, for 
D = 3. This phenomenon, called 'dephasing' in 
is due to the superposition of rapidly oscillatory func- 
tions C(g) on the right hand side of Eq. (p^). In fact, 
for C{q) = ■0fl(a) cos(2u;gt) + 4>i{q) sm{2LUqt), a'^{q) = 
'^q/C'^il) and using y — uiqt — (p- -|- to|. as integration 



variable (m|j 



Am^) one has, for D = 3, 



/3u.(t) 



2^t 



dy CHy) 



X ii^Biy) cos(2y) 



V/(y)sin(2y)) 



(23) 



For large mut the amplitudes iJjr and ijjj vary only slowly 
with y since the leading oscillations have been factorized 
out. This also should hold for C'^{y). Typically these 
quantities are smooth functions oiy/mRt such that their 
y-derivatives are suppressed by l/iriRt. ( For q'^ oo 
we will assume that ipRiq),ipiiQ) and C'^{q) remain fi- 
nite.) The oscillatory behavior of the integrand prevents 
the y-integral from increasing like mut and therefore 
vanishes for t oo. We emphasize that in general the 
asymptotic constant value Am^ depends on the initial 
conditions and is not equal to the temperature depen- 
dent mass renormalization in thermal equilibrium. 
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The time evolution of Am^ for initial conditions 
and a'^iq) = B{q)/G{q), rr? = 1 
and AA'' = 2. B{q) = B is chosen so that the energy ei is 
the same in all three cases. The ultraviolet cutoff was chosen 
to be A = 100 and the total number of modes is 2000. For 
long times Am^ tends to different constant values for initial 
conditions with the same ei but different G{q) and a'^(g). 



FIG. 1. 

with G{q) = l/((7^ + K^m^) 



For a constant Am? the dynamical system simplifies 
considerably and we obtain a free oscillator for each mo- 
mentum g, with frequency Uq . The amplitude of the oscil- 
lations determines in turn the distance to the fixed point 
C* = and — {uj'^)^B^ = 0. Thus for generic initial 
conditions the N --^ oo system does not reach asymp- 
totically the fixed point manifold Eqs. (|l8|). If there is 
a 'dynamical fixed point' for Ato^ the correlation func- 
tions of the system end up oscillating. In addition to 
a{q) the system with constant Am^ has another infinite 
set of conserved quantities, namely A{q) +uj'^B{q). The 
late time behavior of the system is clearly not ergodic. 
Ergodicity is usually thought to be a sufficient, albeit 
not necessary, condition for subsystems of a large closed 
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system to approach the canonical distribution 

The generic osciUatory behavior is also apparent from 
the Hamiltonian structure which governs the dynamical 
evolution of the two point functions of the N oo sys- 
tem. In terms of canonical coordinates Q{q) and P{q), 

Aiq) = a^{q)PHq) + l/Q^q), B{q) = a\q)Q^{q), 



C{q) - ^a\q)P{q)Q{q), G{q) = Q\q) 
the equations of motion become 



(24) 



dtQ{q)=P{q). 
dtP{q) = -ujlQiq) + 
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a^q)Q^q) 

XN r d^p 
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The effective Hamiltonian 



d^q P'iq) 



2 
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(25) 
(26) 

(27) 



2a^iq)Q^iq) 



d^p 



is conserved in time*. The potential is bounded from 
below and the minimum corresponds to the fixed point 
solution Eq. (|l|), i.e. P{q) = 0, u;^Q'*(g) = l/a^(g). 
For any initial value of HcS larger than the minimum, 
the system cannot approach the fixed point and must 
remain oscillatory. 

Let us next turn to the system for finite N and in- 
vestigate the behavior in the vicinity of the fixed points. 
For this purpose we identify u, w, etc., with the couplings 



at zero momentum and therefore 



^0 ' 



B = B{0), 



and so on. In the following we concentrate on classical 
statistics {h — 0). We consider the linear equations for 
small deviations from one of the fixed points (p^), with 
SG{q) — G{q) — G,(g), etc.. For simphcity we choose 
— — 0, whereas B^{q) remains arbitrary such that 
the fixed point does not necessarily coincide with ther- 
mal equilibrium ([l9|). In terms of directions orthogonal 
to the fixed point manifold 



a{q) = 5a^[q) 



]G^{q) 



5G{q) 



1 



\G^{q) 



. 1 
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c.4G(0) 



2 ujIG{Q) -f \So 



(28) 



*This effective Hamiltonian has been found in a different 
approach in Q. We emphasize that in leading order in the 
expansion there is no difference between classical and 
quantum statistical systems and HcS appears in both of them. 



A 



[loIG{Q) + XSoY 



{oj15G{Q) + G{0)5ujI + X5So) 



2^0 



6 



N + 2X 



1 



AT + 8 N + ^uujIG{Q) 
X 1 

1 - 



w 



TV + 8 uu'lGiQ) 
the linear dynamical equations read 



dt5G{q) = -2Giq)c{q), 



dtd{q) 
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dtu 
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dtv = -4 (1 + cr) u + 4uw2G2(0)a(0), 
dtw — SluqU — (3 + 2a) v, 

N + 6 



(29) 
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All coefficients multiplying the linear excitations should 
be evaluated at the fixed point. The excitations SG{q), 
w and s are considered coordinates on the fixed point 
manifold (1^) and do not appear on the right-hand side 
of Eq. (^). ( Note that cJtw^ and dtSo do not involve 
these quantities.) We can therefore discuss the subsys- 
tem {d{q), c{q),u,v,y, z} separately. We first consider a 
scenario, for D > 1, where loq and Sq are almost inde- 
pendent of time due to the averaging out of fluctuations 
with different momenta. In the limit dtUjQ — 0, dtSo — 0, 
the system for x — {a(0), c(0), u, v, y, z} becomes closed, 
dtXi = ^ijXj- On fhe other hand, the anharmonic oscil- 
lator for D = obeys 



dtul = -K^X{N + 2)G(0)c(0) ( 1 
dtSo = -2Ks{N + 8)G^{0)ciO), 



6 



iV + 81 + cr 



(31) 



with = = 1. In order to interpolate between both 
situations we keep k^j and us as free parameters. The 
matrix Aij determines the behavior near the fixed point. 
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For the free system, with A = 0, cr = 0, the six eigen- 
values ^ of A are zL2iujo,zL2iuJo,zLAiu!o, as expected for 
the harmonic oscillations around the fixed point. For the 
interacting system the eigenvalues are determined by a 
cubic equation, 77^ + brj^ + cq + d = 0, in rj = S^'^/Aujq, 
where 
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(32) 



(I + 6S1) 



(33) 



l + CT 



conceivable that the oscillations are damped by the non- 
linear flow - at least for a large class of trajectories. Of 
course, T-invariance implies that for any trajectory with 
flow towards the flxed point there must also exist a tra- 
jectory with flow away from it. Nevertheless, for a non- 
linear flow it seems possible that only a very special class 
of excitations around the flxed point leads to unstable 
behavior. The relative volume of the space of unstable 
excitations ( as compared to the stable excitations) could 
shrink to zero as the number of degrees of freedom goes 
to inflnity. 

We should point out that the general oscillatory be- 
havior found in our linear analysis is not in contradic- 
tion with the damping rates computed in linear response 
theory for classical or quantum fleld theories in thermo- 
dynamic equilibrium H. Indeed, a linear superposition 
of a continuum of frequencies with suitable amplitudes 
can easily lead to a damped behavior, as illustrated by 
the identity J^^ '^zj^^"^^ = e"'^'*'- The damping of 
Am^ shown in Fig. ^ also constitutes a good example. 
It is therefore conceivable that on a microscopic level 
the oscillations are never damped and the system never 
reaches the thermodynamic equilibrium flxed point in a 
strict sense. Nevertheless, macroscopic quantities may 
well approach equilibrium values for i — > 00. 

Let us flrst consider D = where n^j = i^s = P = 
l,r 



jV+2 

2(Af+8)^ 1+a ' 



such that the cubic equation only 
depends on one parameter, a. In an expansion in powers 
of 1/A^ we find in lowest order K = Q, with roots rji = 
-(1 + f ), -(1 + f), -4(1 + f ). This is very similar to 
the oscillations of the free system. However, the relevant 
frequency is shifted to 



The behavior of the linearized system depends on the sign 



of the discriminant, K = ■^b'^c^ — ■§jC'^ 



|5cc?. For > 0, the three roots rji are all real and 
negative, for positive b, c and d. The linearized system 
is oscillatory with eigenvalues — ±2ia;o\/|'7j|. For 
negative K two roots ry are complex and conjugate to 
each other. Since the real part of ^ does not vanish in 
this case one concludes that two eigenvalues ^ must have 
a positive real part. They dominate the evolution for 
t 00. The fixed point is unstable for K < 0. 

The fact that the eigenvalues of Aij occur in pairs 
with opposite sign is a direct consequence of time re- 
flection symmetry of the microscopic equations of mo- 
tion. The correlation functions are eigenstates of time 
reflection T- in our case c,v and y have odd T-parity. 
One infers that the linear system obeys TAT = — A with 
f = diag(-hl,-l,-hl,-l,-l,-M). This implies that for 
any eigenvalue ^ of A there also exists an eigenvalue — f . 
One concludes that an approach to the equilibrium fixed 
point cannot be seen in linear order. For time reflection 
invariant systems the only way an equilibrium fixed point 
could be approached is a purely oscillatory behavior in 
linear order with vanishing real parts of all eigenvalues 
of the 'stability matrix' A. Beyond linear order it is then 



■a/2 = + -NXG, 
' 4 



(34) 



with the temperature dependent 'propagator' {T = B^^) 

-1 



G = 2r ( V m4 + 2NXT + 



(35) 



In order to discuss the stability of the equilibrium fixed 
point we have to proceed to next to leading order 1/A^ 
where 



K = 



N 



cr\3 /47 + 23cr- 13cr2 



1 



o 



1 



(36) 



This is positive for a < cfc 
order relation 



2.982. Using the leading 



(37) 



V2NXT 



one has a < 1 for positive and no unstable mode 
is present. On the other hand, for negative rrv^ there is 
a critical value Vc ~ —4.39, beyond which a exceeds cjc- 
This leads to the most interesting conclusion that for neg- 
ative and XN < 0.026-^, the system is repelled from 
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thermodynamic equilibrium asymptotically, despite the 
fact that the number of coupled degrees of freedom can 
be arbitrarily large {N oo) ! For a given initial average 
energy per mode (~ T) the stability of thermodynamic 
equilibrium requires a minimal coupling strength ! 

In contrast, for = K5 = (I? > 0), we find already 
in lowest order in an 1/N expansion the roots 77 = — (1 + 
cr/2),-l, (-4 + 2cr), with 



K = 



1 



> 0. 



(38) 



The linear system is purely oscillatory. Since X is a con- 
tinuous function of and r (note that r ~ ) 
this remains true for small For example in next to 

leading order we obtain 

+150a + 72) - ^ ^ ^ 

This conclusion also extends to non- vanishing and Kg, 
as long clS f^uj ^ l/N. We conclude that the equihbrium 
fixed point has no unstable direction for linear fluctua- 
tions. By simple continuity considerations this stability 
can be extended to a whole region in the fixed point man- 
ifold around the thermal equilibrium point. 

The discussion of the nonlinear behavior for large 'dis- 
tances' from the fixed point manifold remains outside 
the scope of this note. We only report here briefly some 
results for £) = where Eqs. (0), ( |2^ ) reduce to nine 
coupled ordinary differential equations, which we have 
solved numerically. The stable small oscillations around 
the fixed points near classical thermal equilibrium can 
easily be observed, being determined by three relevant 
frequencies. In this case we have not observed any non- 
linear damping. As the amplitudes of the oscillations 
increase the stable behavior continues for a while. The 
time averaged mean value corresponds now to a 'dynami- 
cal fixed point', which is close to but not part of the fixed 
point manifold. For a Gaussian initial distribution not 
too far from the fixed point manifold ( and for > 0) 
this behavior is approached after a certain time. Finally, 
for large distances from the fixed point manifold the be- 
havior becomes much more complicated. Typically a very 
large range of frequencies becomes important after a cer- 
tain time interval and the problem becomes stiff. Anal- 
ogous behavior has been observed in other anharmonic 
oscillator systems treated in the 1 /N expansion 1 10 1 . We 
will address these interesting phenomena elsewhere. 

In conclusion, the step beyond the leading order 1/7V 
expansion seems crucial for the understanding of non- 
equilibrium field theories. One finds in leading order in- 
finitely many unphysical conserved quantities which pre- 
vent the system from approaching thermal equilibrium. 
They do not survive the inclusion of scattering, beyond 
leading order. New mysteries have appeared, however. 



This concerns, in particular, the existence of a contin- 
uous manifold of stationary distributions within which 
the thermal equilibrium distribution is just a particular 
point. It has also become clear that small deviations from 
the fixed points are not damped in the context of linear 
fluctuation analysis. The understanding of the thermal- 
ization of isolated time reversible systems described by 
flcld theories remains an unsolved challenge. As long as 
such basic questions are not fully understood the use of 
equilibrium field equations in cosmology should be re- 
garded with caution. 
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Forschungsgemeinshaft and the European Science Fun- 
dation. 
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